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ON THE FLOW OF A COMPRESSIBLE FLUID 
BY THE HODOGRAPH METHOD 
II - FUNDAMENTAL SET OF PARTICULAR FLOW SOLUTIONS 
OF THE CHAPLYGIN DIFFERENTIAL EQUATION 


By I. E. Garrick and Carl Kaplan 


SUMMARY 


The differential equation of Chaplygin's jet problem 
is utilized to give a systematic development of particular 
solutions of the hodograph flow equations, which extends 
the treatment of Chaplygin into the supersonic range and 
completes the set of particular solutions. 


The particular solutions serve to place on a rea- 
sonable basis the use of velocity correction formulas for 
tne comparison of incompressible end compressible flows. 
It is shown thet the geometric-mean type of velocity 
Correction formula introduced’ in an earlier paper; part I, 
has Significance as an over-all type of approximation in 
the subsonic range. 


A brief review of general conditions limiting the 
potential flow of an adiabetic compressible fluid is 
given and application is made to the particular solutions, 
yielding conditions for the existence of singular loci 
in the supersonic range. 


The combining of particular solutions in accordance 


with prescribed boundary flow conditions is not treated 
in the present paper. 


INTRODUCTION 


This paper presents a theoretical investigation that 
may be regarded as a continuation of studies initiated 


ms 
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in part I (re pet Ly, in. pare tis sm attempt wes made 
to unify the results of Chaplygin, von Kermen and Tsien, 
Temple and Yeriiood, end Prandtl and Glauert insofar as 
their results were concerned with velocity end pressure 
correction Parone for the correspondence of incompres- 
sible and compressible flows. In addition, two new 
velocity correction formulas were introduced that appeared 
to have a somewhat wider range of applicability than the 
formulas of the efore-mentioned authors. Most of the 
results of part I were obtained with the use of two par- 
ticular solutions of the hodograph equations. These two 
basic solutions correspond to a vortex and a source in a 
compressible fluid. 


It was mentioned in part I that, in order to treat 

the exact boundary problem of uniform flow of a compres- 

Sible fluid past a prescribed body, a general set of 
perticulsr solutions of the hodograph equations had to be 
obtained. Such e study is given in the present paper 
which incidentally helps to clarify the nature of the 
velocity correction factors of part I:= in’ particular, 
the one referred to as the “geometric-mean" type of 
approximation. In addition, many interesting types of 
flows are disclosed from a physical interpretation of 
the particular solutions. A few such solutions have 
already been obtained and discussed by Ringleb (refer- 
ence 2 


Several mathemetical approsches exist by means of 
which particular integrals of the hodograph equations 
may be obtained. Two such approaches, mentioned in 
part I, may be attributed to Chaplygin (reference 4) and 
Bers and Gelbert (reference ) and are snalogous to an 
exponential and to a power-series approach, respectively. 
Another method of defining particular integrals is the 
integral-operator method of Bergman (reference 5). In 
the present paper the differential equation, first used 
by Chaplygin in his treatment of jets (reference 4), 

provides the besis for the definition of « complete set 
of perticular solutions. 


The stope of the present paper is limited chiefly 
to a systematic study of the fundamental solutions and 
to the physical interpretation of some of the particular 
flows represente d by them. The combining of particular 
solutions to represent uniform flow past e prescribed 
body is not treeted herein. It is believed, however, 
thet the present study may serve as a basis for further 
development end clarification of this importent problem. 


NACA ARR No. L129 3 


a 


99 


For 


SYMBOLS 


rectangular coordinates in plane of flow 
magnitude of fluid velocity 


angle included by velocity vector and posi- 
tive direction of x-axis 


eensivy of fluid 

pressure in fluid 

weloecity of sound in fluid 

Mach number (q/a) 

quantities referred to stagnetion point q=9 
velocity potential 

stream function 


ratio of specific heats (approx. 1.4 for air) 


(approx. 5/2 for sir) 


maximum fluid velocity (corresponding to 
B= pS = 0) 


dimensionless speed variable 


ai q ie Mi 
7 Ss TS ees 
2 ves 
( 2Ba, 2B + M 
; i we us ; 
Songe Value or 7 tj = = Bp Prexs 


ee DN, ame 8 
1/6 for eir) 


(or 23 ee Sy, Phe anes Of Tia 
Aan 
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GENERAL PARTICULAR SOLUTIONS OF THE HODOGRAPH EQUATIONS 


The linear equations in the hodogrsph veriables @ 
and q, which relate the velocity potential 9 and the 
stream function w for the steady two-dimensional flow 
of a nonviscous compressible fluid, are 


“ Ov 

mpsS. 
i 

a 
a 


¥ 
<< 


> (1) 


ote 
al 
tl 
\ 
> 
N 
Ke 


r the adiebatic equation cf stete between 
a 


7 4. 
Low 
Luvs 


end 
(See 
In the incompressible case. T-—90, equations (1) 
can be expressed in the Cauchy-Riemann form. ~Perticuler 
solutions "2 = % + iW can be expressed in this case es 
any anélytic function of the complex variable 
Wie o-$;4 ters (2) 
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or as any analytic function of the related exponential 
function 


FAN = qo? 9 (3) 


Thus, an infinite set of perticuler integrals of equa- 
tions (1), in the incompressible case, referred to herein 
as “the powers set," is wt. When k is a positive 
integer, the particular solutions vanish at the origin 

(@ = 0, log q = 0) and, when k .is a negative integer, 
the varticular solutions are infinite at the origin. In 
the case of nonintegral values of k, the origin is a 
branch point of the functions w*, 


e€ 


Another infinite set of particular integrals of 
equations (1) in the incompressible case, referred to 
herein as "the exponentiel set," is 


(e7iw =, gen ihe 


where, again, k can take on any value - integral, non- 
integral, positive, or negative. 


In the compressible case, the particular solutions 
Corresponding, to the powers: set wk _ (that LS ,.4 the pars 
ticular solutions which reduce to w+ in the incompres- 
sible case T->0) depend on whether the coefficient of 
w is real or imaginery - a consequence of the fact 
thet, in the compressible cese, % and w do not 
satisfy the same differential equation. For example, 
for k=1, the two functions corresponding to w 
and iw, which have been developed in part I, are 


W= 868+ iL 
and 
iW = 1(6 + iL 
where 
b= toe gq + f(t) 
and 
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and :f(t)--and £479) cceach) van ieh) fom + 4S) Obsd See 
equetions (26) and (27) of reference 1.) 


The development of other functions corresponding to 
power set wk, for positive integral values of k, 
ws according to the method of Bers and Gelbart. 
expression (22) of reference 1.) Since the present 
efly concerned with the functions corre- 

the exponential set e~t4W, the powers set 
ier discussed. 


Cheplygin Differential Equation 


Th +4 3 iy 
The Functions P, and Qy 


Corresponding to the exponential sets in the incom- 
pressible case 


are Sm: 
2 q Sim ke +e “eos ica 


there appear in the compressible case functions designated, 
respectively, 


Py (q) cos k6 - 1 Q@(q) sin ke 
and 
Py.(q) sin k@ + i @(a) cos ké 


where the func q) “and Q@, (a) satisfy second- 


P,,( 
21 equations. These equations are easily 
ituting in equations (1) the product- 


Fa) 


+ 

u 
order different 
obtained by s 
type solutions 


~. 


f= Pla) Say, (ee) 


¥ = @ ta) S20 (Ke) 
< Ki cos 
k Kw ; 
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In view of equations 


(1) it is observed that 


O dQ, ( ) 
ky P.( 4): = - 
k 0 aq 
(5) 
dP},.(q) a {P 
ey er ie 
aq dq eg 156 
The functions @, (4) satisfy tne second-order differ- 
ential equation ~ 
PG Se. AQF) ess 0. 16) 
M JQ dae -~ kK \i- M Qk = 
The functions P}.(q) cen be obtained from @,(q) by 
means of the first of equations (5). Equation (6) may 
Berrecucead to 2 standard type by introducing + as the 
independent veriable. Put 
ata) = 4° ¥, (47) (7) 
a4 qQ k 
where clearly Y, Ae )—>1 es t— 0 (incompressible 


case). 


With ne “use of 


the symbolic relations 


cle a 
(Seen ice 
7 aq at 
2] Ps 
ae sheet oa oF 
7 ee 
a age = arc i at 
and the relation 
uz 2 28T 
ven, 


the desired differential equestion is 


ay, 


T(L-7T) ——— 
aT 


+ [(c+1) - (e+ 1 8) 


d¥. 
t| = ae +e pk(k+1) ¥,=0 (8) 
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Equation (8), which is of the hypergeometric type, was 
first introduced by Cheplygin in his memoir on gas jets 
(reference 3). 


The Funetions %. “ahd Yop 


Chanlygin treated the subsonic flow of a compres- 
le fluid through jets with straight-line boundaries. 
such problems the hodograph variables § and q 
natural veriables in the sense that the solid and 
id boundaries ere described by 8 = Constant and 
Constant, respectively, and only the particular solu- 
of equetion (8) with positive characteristic 


oct 
2) = 
ct eae ©: Yl 2 @ ty oO 
ei 
Sd 


ndex k are needed. \In the present, paper ‘a compleuc 
ordered set of particular solutions of equetion (8) ts 
obtained, which extends the results of Chaplygin into 


a 
the supersonic renge and to negative values of the 
index k. Two types of solutions of equation (8) for 
nonintegral values of k are 


Wilt) = F Gye. Pyesektlsag) (9) 
and 
Tut). s TK PQ -k, by-k, 1-k; 7) (10) 
where 
Qe aug do, = ie Ta 6 
k 
ayby = -5 (ke + 1)8 
and 
} ¢ oh ip Vase 
F(a, by ¢; 7) & 2+ S293 ane )b { dey oO 
Cc 2 Ue ea 
It is now shown that only one of the solutions need be 
used. For positive values of k, the ea es that 
Y,, (0) = 1 excludes the use of equation (10). For nega- 
tive values of the index, the solution @_(a)=q7 Yo, tt) 


obtained with the aid of equation (10) is, except for a 
constent fector, equivalent to the solution 9. (apsqk WAF) 
obtained with the eid of equation (9). Thus 
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foiay sa” T., (4) 


q7krk F (atk, biytk, wih; 1) 


| 
4Q 
a 
4 
ye 
ty 
fas 
ity) 
aa 
ion 
a 
Ne 
oa 
+ 
ht 
4 
be 


| 
~~ 
WD 
o i 
12) 
ine) 
e 
a 
-D) 
s 
~Q 


Hence, only the solutions given by equation (9) are 
needed for the determinetion of Qy- (4) and Q1.(4) 
Then : k 
Q,(4) = gk Hyatt) 

= gk P(e Dy, fkt1; T) . Fany 
end 

-k 

mtakee.d ¥_, tr) 


Sia" F(a), by, -ktls T) 
= qk P(a,-k, by-k, -k+1 1) (12) 


Cbserve that both types or hypergeometric functions 
appearing in equations (9) and (10) are utilized in the 
expressions for @ (qa) and Qj(q)- 


The foregoing discussion has been limited to non- 
integral values of the index, positive or negative. When 
the index is integrsl and positive, equetions (9) and (11) 
remain valid. When the index is integral send negative, 
however, equetion (12) does not in general lead to a 
meaningful solution and consequently another independent 
solutsom @is+ to be sought. The desired solution fer 
Y¥_1.(1T) -in»such cases contains ea logarithmic term and 
again is subject to the condition that it reduce to unity 
for T= 0 (incompressible case). The expression for 
Q;,(4) is then given by 


CC) ae nae ae (13) 
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Coma 


*so[duexe TeTosds Zurmop~{[oy e9 Jo swos UT pszBAQSNTTT St Fusuaqguys stu, *z[nses 

a4Tuly pus ZUsAeTea wv pretA you seop (21) uotzenbs gt A[UO paezTTyqn eq oF St (Lt) uotgenbe 
qeuy ‘toaomoy ‘yno pasqutod eq prlnoys 41 *TetwmouxkTod 8 saulioceq eine? uoTZOUNF 944 puR O19Z 
stenbes © 4uvgsuoo suzy ‘*H *°* “Z *‘[ sontea oy jo Aue uo sexey Aq “ag %y Jt ‘4849 e20N 


PCL = 45) 9s yes 
Gam) GG IO) GE) 
pus 
er ae. Oop ee - ee @ + a)(1 + %) t2 
ae A - -=---- -_+— + = +— : ah 
Z+4 t T et L ie L I (1 . iq) iq (1 O Hele 
= = 
a ioe tt ET g0t (x ete 7 ‘te)d Ljo + 
T+ it Hoek oe qh . A 
PG. =)-f = 3) 
NS naw? $$. |, (1-) + 
a Ce ree CI) (CID Ce CREED (CD 
(¢-a)(2 -a)(l - 4) Pe 
Datel ee li eee : os 2 hk nn sree mS 
BA AB a a a )(t + i= iq) y - a )(z +H- ta )(t +i = %e) (4 - %e) 
bs : - 4 os ~ 5) ——— oe a Bs Tine: 
(t +4 - Wa - a(t +a - 7a - =) (4 - *a)(a - 78) eas 
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Cease of y= -l: 


nq Sea ¢ r . 
Consider as an example the von Karman-Tsien treat- 
ment of compressible flow (reference 5) in which the 


adiabatic index 4 = -l or 6 = 5. Then 


k+l 
ak 
by, = 5 


For a negative integral index, equation (13) may appear 

to be appliceble, in which case the expression for 

Pont) would be a polynomial of degree k= i. An 
exemination of equation (12) shows, however, that for this 
ease no infinities arise and that, when the index is 
negative, integral or nonintegral, 


et = r(4sk, -§, lek; ) 


The hypergeometric series represented by Yo, con- 
verges for velues * 0 S lt| < Ll. For the present case of 
w= et or p= os values of T corresponding to 


positive values of M 


aE outside the range of conver- 
Bence. A closed express 

at) 

FS 


Q 

eon for Y¥_4(1) ean be found, 
ich, by analytic continuation, 

it wetlges*or Ard PTs 


r 


however, for this cas 
iS uneretore valid for 


— nu l-k ee Fa ) 
Lar ai ) = FP sk, 2 ’ 1l=k 3 T 
C 5 | kk 


k #j(1 = nee 
a 


Similsrly, from equation (9), when the index is positive, 
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1 


2 
Observe that 
Al 
Q,(4) = 

Z Q_).(a) 
k 
: | 

= |9————_——— (1).) 

Lhebi neta oul 


1 ee ° 3 4 4 a 
This identity for the von Karmen-Tsien case corresponds 


to the identity qk = + for the incompressible case. 


q 


Gage on te ve 


Then, for the positive index, 


Q, (9) =.9,%,.004 
ic) BCL, -6, 28.0) 
J SRE at 
= Q Eo Ae Ces, {255 
CB ts 


For the negative integral index, it may appeer at first 
glance that equation (13) is needed; however, equation (12) 
does yleld a relevant and finite result and accordingly 

is the equation to be used. Thus 
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2 2 
Lim F a,~k, by-k, 1-k; T) = ks EB T - Br ae, 12 + PUB) 23). oo 
k>1 e ox ay 2x3! 


and therefore 


Q_4 (7) sat {ied eBettha a. ne] (16) 


wasenoL. Kk = .O:: 


MWhe exceptional ease of k= 0 is directly treated 
by means of equation (&). The differential equation for 
YotT) ap Got). then Is 


he general solution of this equation can be written as 


HH 


T _, 
5 3 aT 
Q, (4) = Cy log q + Cy, i [a - rf - L| ad Co 
Vo 
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where CC, and Cy are arbitrary constents of integra- 
tion. The constents C, and Cy are determined by the 
imposed condition that the expression for Qo (a) reduce 
in the incompressible case simply to log q. Then 


xd 
a=" Ss 
Co = 0 
and therefore 
a: i He velty 
ala)=togatz f fa-r)P-a} 2% aq) 
JO 


In a Similar manner, from the differential equation for 
oa 
fe) 


a | 7h a tP Ary 
dv | bees We ee an 


the expression for P, is obtained as 


Te 
Pika? =.log q * 2 | eer = | _ (18) 
-T 


It is remarked thet the functions Q,(q) and P,(q) 
are identical with the elementary functions L(q) 

and L(q), respectively, introduced in part I (refer- 
ence 1) and are associated with a vortex and a source 
type of flow. 


The Functions Ry and Si 


A linear homogeneous differential equation of 
order n” can, in general, be reduced to a difficrential 
equstion of order n- 1 by means of an exponential- 
type substitution for the dependent variable. Chaplygin 
mede use of such a substitution to reduce the second- 
order differential equations satisfied by Py. and Q, 
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to first-order equations of the Riccati form, in order 

to study properties of the functions PF, atid Q.° in the 
subsonic range for only positive values of k. In the 
present analysis the Ricceti equations are also found 
useful in order to extend the study of the functions Py, 
and Q to the supersonic range for both positive and 
negative values of the index -k. 


The second-order differential equations for Py; 
end Q,, with t+ as the independent variable, are 


tee sae eo Me “ae 
tee ee oF 3) Tt = 


The corresponding first-order Riccati equations are 
obteined by substituting for P, and Q, new dependent 
variables Rk, and S,, respectively, es follows: 


k 
Lee, aT 
Pose 
or 
Sea 
sf a ie 
et |.d : 
z sla 1B me (19) 
4m u 
and 


| 
@ 


(a OS 
he 


16 NACA ARR No. L4129 


or 


| 
= |Y 
. t= 
Q. 
PS 


w | 
ior 


—— Loe e (20) 


The equations satisfied by R,(t) and S;,(r) are 


ate 1+ (2psi)r 6 R Sb ieee aearseae (21) 


E RY Se 2 
dr L=t{2p tt leit QT Loses. 
and 
dS, 6 k 2. De ee eee 
ae ae eee 2 Ce ee eee oe 22 
at te oe em. s, Lae as 


Initial conditions for R,(T) end S(t) are found by 
examination of the incompressible case t-—~>O0O. In this 
k 


cese Pr = Q@ = 4 and, Silmee . et re = -q an it follows 
ay q 
from equations (19) and (20) that 
Ry, (9) = S;-(0)° = 1 


Jpal 


The following important relation exists between the 
functions Ry,(t) and S(T): 


Ry, (7) S08 b= 


| 
ke 

i 
is 


I (23) 


Equetion (23) cen be verified directly from the hodograph 
equations (1). It may be noted et this point that this 
result ispof significance in connection with the 
geometric=mean type of velocity correction factor intro- 
duced in pert I and is discussed more fully in a later 
section. 
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Before the functions R,(t) and S,(T) are treated, 
certain generel observations can be made regarding the 
functions Py(r), @,(T); R(T), end S,(1T). Chaplygin, 
who limited his investigations to the subsonic renge and 
to positive values of the index k, hes shown thet a, 
end consequently the other functions possess no roots 
for any value of the independent variable in the subsonic 
range, with M= 0 excluded. Inthe supersonic range 
a P,.(T) and Q(T) in general possess zeros. 
Certain reletions obtained by means of equations (19), 
(20), and (23) between P,, Qo. Ry, end S, at the 
zeros of P, and Q, are summarized es follows: 


1) | Mex or min 


Vics or min 0) 


It is remarked thet the number of zeros of Q,, 
as a function of the index k, can be found from an 
expression developed by Klein and Hurwitz (reference 7) 
in connection with the zeros of the hypergeometric func- 
tion. In generel, the number of zeros increases with 
the pagent tude Ofimtaeé timdex tk cand*is Infinite efor 
k= 2a, 


A further observation of interest can be made in con- 
nection with equetion (24). Cheplygin has shown that, 
for positive finite velues of k (and the same is true 
for negative finite values of k), the functions Si (T) 
are not zero for the sonic value T=T, or M=1. 
From equation (23) then, it follows that the functions 
R(T) =O norm NM = a. 


In view of the relation between the functions hy, 
and §, given by equation (23), only S, need be dis- 
cussed. The Riccati equation (22) may be used to discuss 
certain properties of the function S, but in general, 
for numerical evaluetion, the original definition (equa- 
tion (20)) in terms of the function Q, may be used 
girectly 3 


18 NACA ARR No. LyI29 


, 26 1. oo 
k k Q, aT 
or 
dy; 
ge 1 theo 
Yr at 


In general, the functions S, are expressible in infi- 
nite series. For several values of k, however, 5S, 
can be expressed in closed forms. For k=O and 

k = fo, S; may be obtained by a limiting process from 
equation (20); however, for these special cases the 
Riccati equation (equation (22)) yields the results 
directly. Thus 


So = (1 = T)P (2h) 
and ) 
eee (26 a9 ive 

+00 wae" 
ve 
eitie Se" (25) 
he cases k=1 and k= -1 may also be ex#ressed in 


closed form. With the aid of the equations (15) and (16) 
for Q, end Qj, equation (20) yields 


EIS (ha Bet =r) 


S.. = KR wi2 (26) 
: Lb =O (2 -7)Ptt 
and 
ee he (27) 
a ey ae 1h] 
Be ie 


In order to illustrate the behavior of some of the 
functions thus far introduced, a number of tables and 
figures are given. All the calculations heve been 


— 
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performed with the adiabatic index » = 1.l. Table 1 
gives values of Y, as a function o Mor “gyetor 


several positive and negative values of the index k. 
Figure 1 shows the Y, functions plotted against M. 
Values of the 5, and R, functions are given in 
tables 2 and 4 and are plotted against M in figures 2 
and 5 ° 


The Functions f,(T) and 2,47) 


In the incompressible case, the sets of functions Q, 
Pnas) fj, Can be reduced to a single function log q by 


r ae 
means of a simple operetor Z 1° Thus 


and 


This seme operation applied to the functions Q, and Py, 
in the compressible cease serves..to..define.twousefulb sets 
Orvmemctions iog q+ f(T), and, log.q,t g(t), respec= 


tively. Thus 


Oe ie : 

pers &, = logy *of),(q) (28) 
and 

1 Sip we i deat “e.(er.) (29) 

a E k ae 
From equation (7), namely, 

1 
it follows that 
f(t) = = log ¥y,(1) (30) 


z 
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From equation (5) for Py, and equetion (20), which 
defines S,, 


Al; 
P = — 0,5 
k Pia 1)P k k 


It follows that 


1 


ry |r 


log Y, (7) S)(7) 


£:,( 7) 
i (1-7) 


Si 


f,(7) + = log (31) 


(1 - 1)? 


For example, for k = 1 and k = =] and with the use of 
equations (15) and (16), 


De Fess (Gtd Yom 7) Stl 


£3 (7) = ‘loge 
(Bn Lot 
pie ete Be =i 
g(r) = log Ga eat (33) 
(Bo EGC = iF) 


2-6-4. - 


f_4(me- 10¢ 43 = ee [2 -(1- aa (3h) 


(38 +1) - B+ Tt) - 7)P 
2(p + 1)(1 - 7)P 


fis 


.i{T) =- log (35) 


For k=0 and k = tm, equations (40) and (31) require 
a limiting process for their evaluation. Alternate forms 
for f,(tT) and g,(T) may be obtained, however, by 
means of equations (19) and (20) defining Ry, (7) and 
S(T), which yield the results. for k= 0 and +k = a 
dipsctly. Thus 
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Ste jah ee 

ry Pe ee | at 

P(t) 25 | [sy.() 1 {St (36) 
WO 

and 

OT “6 

a(t) = 5 | (R(t) - 1] (37) 
WO 


where R,(T) and S,(T) ‘are related according to equa- 
e207 (25). Then 


yr =] 
i H aT 
f(t) = 5 [ca auadt j= L poten (36) 
pit 
rr 1. - bee 
is H US roe hate Fy) Var 
e(t)=$} | toes lN je (39) 
Pa A ee Oke | in 
pe a «tte t) | 
and 
Or (— Hiv? j 
, Pe eS eae aT 
2.785) Sr, = f= L eye 
toT) = GaorT) an af Pega : (10) 


it is worthy of special notice thet the functions f(T), 
eat i> and Eben (T) are identical with the functions 
f(t), -el(t)s. aad hi*«), respectively, which formed the 
bas2s,of part: I feed gwaree Lj» »Iin additiony, the i el 2 
Siehsypdegeast £,{1)5  ,legia + Bit), _end legndg hips (T ) 
are identacal with GhetfunetLrons® ££, bij. land: 1 Hy 
respectively, which weie introduced in pers I. 


A numoer of functions f;, and Ex have been 
calculated, with = 1.4, for several positive and 
negative values of the index k, and the values are 
given in tables | and 5 and plotted in figures ) and 5. 
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The opportunity is teken here to note thet, for the 
von Kerman-Tsien case (y =-l or B=-= 


i 
t 
hb 
ie) 
09 


fy = && 


1! 
a) 
e) 
t0)¢) 
ine) 
= a 
i] 
a 


and that the sets of functions P, and Q,, as in the 
incompressible case, are reduced to a single function by 


the operator = log; namely (compere equation (1l)), 
avi = Ne 


Lut yh = yA 


In fact, the complex flow potentiel 9 + iv can be 
expressed as an anelytic function of a single complex 


aJ/L - Me 


veriable § + i log q ——~—————. Tsien has made use 
1+ - Me 
of this complex variable in his hodograph treatment of 


he compressible flow past an elliptic cylinder (refer- 
ence 8). 


log q + Log 


Velocity Correction Factor 


The solution of the problem of an exact corre= 
spondence between the flow past a prescribed body in an 
incompressible fluid end the flow pest the same body in 
a compressible+fluid is a: Giliiteult matters (rhss 
problem cen be solved exactly for certain types of flow 
patterns (not pest closed shapes), such es flows inside 
or outside angles or channels, and for certain flow 
Singularities such as a vortex, source, and doublet = 
types of flow which cen be associated with the particular 
solutions Q,.- Some of these types of flow are illus- 


trated by examples in the following section. Combining 
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particular solutions to represent uniform flow past a 
prescribed body is e complicated process, since the treat- 
ment of infinite series in the functions Q, for both 
positive and nezetive velues of k is involved. Further- 
more, the process of returning to the physical-plane 
eriables from the hodograph-plane variables hinges on 


nonelementary parts of differential geometry. Certain 
types of jet at eee can be properly treated in the sub= 
Sonie renge by series in ¢ with k positive, as was 


shown by Chaplygin (reference 3). Thus, it appeers that 
much work remains to be done in order to render feasible 

xact and practical solutions for uniform flow past pre- 
Scribed bodies in a compressible fluid. Beceuse of the 
difficulty end complexity of the general problem of flow 
in a compressible fluid, attempts have been made by a 
number of investigators to obtain results by means of 
Wetocity correction formulas that serve to place in 
eorrespondence velocities in an incompressible and in a 
compressible fluid. 


In part I the velocity correction factor was dis- 
cussed c perticular reference to the two functions L 
end L and P, of the present paper r) associated 
with a a ae and source type of flow, respeciively. The 
main justification for the results of part I was the 
yielding and the unifying of the results of Chaplygin, 
von Karman and Tsien, Temple and Yarwood, and Prandtl 
and Glauvert. The knowledge of the infinite set of func- 
tions Py, and @ discussed in the present paper can 
now serve to Laces further on a reasonable basis the 
concept of a velocity correction formula. 


In order thst a single velocity correction factor 
be feasible, even for a flow associated with a particular 
solution, it. is necessary that Py =Q,» Consider, for 
example, the functions Q, and P, insofar as the first 
power of the veriable tT is concerned. It can be shown 
easily that 


yy fr 
fe 
(e) 
09 
B) 
wf 
I 


log q + f(r) 


v 

bh 

(e) 
03 
Ke) 

1 
hoje 
i?) 
+ 


2); NACA ARR No. Ly 129 


5 
it 


a log Pie = log q + g(t) 
jy 
= 1 a 
oe a - 5 Bt 


Thus, to the first power of T and independent of k, 


fplr WS appe) 


Ww 


ue 
“3 PT 


Py = Q 


2} 


4 


The neture of the correspondence between the incompressible 
flow and the compressible flow is such that 


rv 


Saath ah 
Va Fee 


Without going into any details here of the field point 
correspondence or of the boundary distortion, the veloci- 
ties in the incompressible end compressible cases may be 
placed in correspondence as follows: 


: aap 1 
(log q)y = (Log d= Seay) 


c 


GQ; = ae (112 ) 
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This result implies thet the complex variable 

; Ht : : : 
@ + i (log q - 5 BT in the compressible case corresponds 
to the complex varisble 9+ i log q in the incompres- 


sible cese. Equetion (2) represents the approximation 
ef Temple and Yarwood discussed in part, I. 


Consider now the functions Q, and ie inscfar as 
iareo velues of the index k are concerned. It is 
recalled thet, as the index k —>to, 


and that 


where 


Then, es k —>to, 


rd 

ws 
u 
£ 


ae 


The function h(t) is expressed in integral form in 
equation (1,0) and hes been eveluated and tabuleted in 
part I. (See also table lh and fig. of the present 
paper.) The correspondence of velocities in the incom- 
pressible and the compressible case 1s given by 


R 


qa, = qyet) (43) 
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Equation (43) constitutes the geometric-mean velocity 
correction formula introduced in part I and is limited 
to the subsonic range 0 S<m$1. It is observed that, 
for positive values of k, h(t) lies between fy (7) 


and g;,(T) in magnitude. Moreover, the deviation of 
£44 g,( 1) 
from e xf) and e < 


(T) 
ubsonic range. (See table 6.) 


is quite small in the 


The foregoing remarks, together with the fact that 
the geometric-mean type of approximation contains the 
results of Chaplygin, von Kdrmdn and Tsien, Temple and 
Yarwood, and in the limiting case of small disturbances 
to the main flow the exact Prendtl-Glauert rule, lead 
to the suggestion thet it may be adopted as an over-all 
type of approximation in the. subsonic range. 


Flow Patterns Corresponding to the Particular Solutions 


Before the flow patterns corresponding to the par- 
ticular solutions pi. end WY, given by equations (hh) 
for the compressible case are discussed, it is instruc- 
tive to examine the incompressible case. Consider the 
complex velocity potential 


Q 


g+ iv 


eae (ih) 


where U and n are constants and 2 = x + iy. itis 
well known thet, if n= tT where a is an engle between 


Oe: 
O and 27, equation (i) represents the flow in a sherp 
angle. For example, the flow inside a right angle is 
obtsined with n= 2 and the flow outside a right angle 


is obtained with n= ae Again, the value n=] or 


or @= etm corresponds to the flow around a semi- 
infinite line. Clearly, all the gle flows are obtaine 
with values of n between 1/2 and oo. Other types of 
flows are given by other values of n. For example, 

n= -l corresponds to a doublet and the remaining 


2. a > 1 1 
a= T corresponds to a uniform flow and the value n= > 
d 
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negative integers ere associated with singularities of 
higher order than the doublet. In addition to the flows 
described by the powers 2", there are the two funda- 
mentel flows, the source and the vortex, associated with 
the function log =z. If, now, it is desired to obtain 
generalizations for the compressible case of the fore- 
going perticular flows, the procedure is first to express 

or wW for the incompressible flow as a function of 
the hodograph variebles q and 6 ana then to replace 
a. Oy Py. or GQ.» respectively. Several examples will 
best illustrete this precedure: 


(1) Consider the compressible generalization of the 
angle flows. By meens of the relation 


an <4 
as 6 TW 


where w= 8+ i log q, the hodograph complex variable w 
is introduced as independent. variable, in place of 4a. 

o \ 
From equation (ji) 


60 42 2. ns 
— = nUz 
Big 
=) iw 
Hence 
de 
+o \re L 
lee us (--i”) 
if 
ne 
(Un) 
and 
SS 
Vim 
T -i1W N= 
a= —2— (27) 
n-1 


ia) 
(ee) 
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Then 
eds 
oF U n=l a; 
ee (Serer Se a Sin =~ 8 
n-1 
(Un) 
Lf one is replaced by k, the compressible generali- 


zetion of the angle flows is given by 


J 


Wy = - ha tad Q, sin ké (45) 


The inside angle flows are given by values of k in the 
range 1<k<® and the outside angle flows, by values 
of S4e 9 in “thie? rence a < =k < @, For example, k=2 
for the flow inside a right angle, and k = -2 for the 
flow outside a right angle. Other types of flow are 
Given by values of k in the range -l1< kil. 


The case k=1 or” n= im: is exceptional and, in 
fact, corresponds to the incompressible flow 
os es 
Q =e (46) 
where coc is a constant. 
(2) Consider the compressible generalization of the 


doublet. The complex velocity potential for the incom- 
pressible doublet at the origin is 


‘o) 
iH 
NI 


The reflected-velocity vector is 


a0 2 on 
az Zo 
-iw 
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Hence 


ana 


whe streem function for. the incompressible doublet is 
then given by 


The compressible generalization of the doublet is 
therefore 


(4) Consider the compressible generalizetion of the 
Source. The complex velocity potential for a unit source 
at the origin is 


The reflected-velocity vector is 


dO wl: 
az «2 
om e7 iW 
Hence 
zZ = 6lW 
and 
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ty potential for the incompressible source is 


3 
3) 
(o) 
< 
(49) 
aa) 
1o) 
(e) 
= 


¢=-logq 


compressible generalization of the source is then 


(l.) Consider the compressible generalization of a 
int vortex. The complex velocity potential for a 
vortex of unit strength at the origin is 


isos! Ploeg 


The reflected-velocity vector is 


aot 
dz Z 
=) 67 tw. 
Tience 
z = =1el¥ 


end, except for an additive constant, 
Q=w 
The stream function for the incompressible vortex is 
w= log q 


The compressible generalization of the vortex is then 
given by 
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Trensfrormation from the Hodograph 
to the Physical Variables 
Given the velocity potential @ and the streem 
function wW in terms of the hodogreph veriables 6 
and qd, it is possible to express the coordinates x 
end y of the physical plsne in terms of" and ‘q. 


From the basic flow equations 


ge! oy 
og _ _ Po oy 
oy p Cx 


it follows (see equation (5) of reference 1) that 


cs aN 
ax = i (2 cos 6 - ake ea J sin 8) dq 
oe | 


vi Ce sin 9 + Sraae cos ¢) an) 


Equations (4.7) relate the differential line elements in 
the physical x, y plane and the hodograph 8, q plane. 
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When expressions for ZY and w as functions of 9 

and q are known for a given flow, the integrals of 
equations (7) are the equations of transformation of 

the: @y 0d) tebe x, Jy coordinates. It may be remarked 
at the hodograph flow equations (1) are the integra- 
bility conditions for the differential equations (7). 

The right-hand sides of equations (l17) are therefore 
perfect differentials. 


Consider one set of particular solutions from equa- 


tions (i) 
% = Py(q) cos ké 
wW = - Q,(q) sin k®é 
where k= 


Ll and k = 0" “are 6xeluded. By the use ef 
); 


t 
(5 it can easily be verified thet 


equations 


oe Po g \ cos(k + 1)8 
q rene! k | a ae 


a 
2 
7 


P : eee 
: ( Y ®2 a, odd oe denis Constant 


(48) 


k Py, 0 sin(k + 1)8 
ee <2 ee 
"kk alc 2 0) k+l 


P fe) Sin(k = 1)9 
Ss k ops pel fala toa iat Jad 
& + ae ae | + Constant 


The equations of transformation corresponding to the 
other set of particular solutions from equetions (l) 
obtained by replacing in equations (48) the cosine by 
the sine and the sine by the negative cosine. 


The excluded ceses k= 0 and k= %1 are now 
trested. For k = 0, one set of particular solutions 
corresponds te a source and is 
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2. ae Ee 


f= 6 


Equations (li7) then yield 


x = — cos 80 
oa 
Po 
y= — sin 8 
* og 


ghe other set of particuler solutions corresponds to a 
vortex and is 


Equations (1:7) then yield 


Pa 

il 
RO} Te 

07) 

bre 

D 

ram) 


home kes lk with 
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sity Aa e 
ee L ad es Po 2 Sar eG, ee ‘a nt + Po@ 
Sa ENG erie e\4 pa 


With the use of equations (15) and (5), 


ae ‘AS (dee 
Ky = 5 1 - cos 26.4 les gq 
(6 + 1)r(1 - 7)? 


| one 
{1 
Pope 


where g(T) = B(7) by equation (39) and is evaluated 
in part I (reference 1). 


For k= — i) witch 
Py = Pl Sin 6 
Wy = Q cos 8 
4 set nies Brl 
i =|) Eee sin 26 + 6 


Fie = Vee StL ae eee, eos 28 + Loe >q (50) 
(6 + 2) eC = 
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For k= -l1 with 

#_1 = P_j cos 6 

vy. = a4 sin 6 


equations (/17) yield 


a! 
eS Pe » Po*s)} og + 2b | y OP e “1 \ 
X_y }cos 2y q 
> A ea / 2.9, 68 pe cay 


With the use of equations (16) and (5), 


x a | af + 2 pe Ae oe ae (1+ 67) | ¢oa 2a 
eS hg R+1 (1 LryP Gd 
= Ae 
eee ioe qe 2 = g(r) 
hoy : 
ha, (6 + l)a, 
ae ee hs - 
eee ge |i = TP (51) 
Tyas OPO Sy ad et sented. + pall sane Bo 
et (t= 4)” p+ 1 
i 
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\N 


Hom) densa wavtin 


ay = Fy an 8 


Vy = Qj cos 8 


$7 3s 26+ 2 aig Sa 8 (lL + BP) sin 28 
Wor if pel (aera Bt 2 
Wipe aay 
ha © 
wie: 
ee eee weeny MOE ree bys Wester: (52) 
-1 Wee nae . 6 2 
HQ f + (1 T ) pwd 


+ he log q + — ee ater) 


ha.“ KiB + Le, 


4 eee SE et i - 1 
8(p6 + 1) ba,“ (ie 1)° ‘ : 


Ringleb (reference 2) gives an example of the flow 
of a compressible fluid around a semi-infinite line. An 


‘ ; iy Nie oe : Beh a 
examination of Ringlebts stream function v= q sin 98 
shows that it is a linear combination of W, and W_43 


that is, 


2 


ph ad OR $ 
v= é eae Qh + a1) sin 9 
~ nrc # 
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In fact, all the externel angle flows Gi< -k < o ) are 
nonunique; for, in view of the discussion preceding 
equation (11), a general form or W_, is 


J na ike i tee ‘ 
ee =q E Ye) + v(7) | sin ké 


where A is an arbitrary constant. 


Observations on Limit Lines 
In the Bee. section there are reviewed briefly 
certain conditions, discussed by Tollmien (reference 9) 
and Ringleb (reference 10), with regard to possible limi- 
tations on the potential fiow or en adiabatic compres- 
Sible fluid. 


Consider the family of streamlines 


W(@2q) = Constant 


Then along a streamline 


and, from equations (1:7), the line elements along a 
streamline are 


Be 
ah | 
ex = 22 [200 - 1)(2)*- "| Bb ae | 
sf ‘ \ (53) 
E 7 
ay =o | 20 - Gy - GO" ES as | 
f°) — a v/a 


Singuler points elong a streamline are cheracterized by 
the vanishing of the common factor of equations (53): 
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(quy° - 3 - 1) = (5k) 


(Stagnation points at which = and : vanish, the 
q 
vortex for which et = 0 and the source for. which oy = 0, 
g oq 


are excluded from this discussion.) Observe now from 

equations (7) that the Jacobian of the transformation 

from the hodogreph variables @ and q to the physical= 
plane veriables x and y is given by 


GD = 23h u- ew) 


e54q 


- “zur -1 cn (55) 


1 
foul 
> 
ob: 

ft 
Sys: 
Cy 
ele 


Thus, the vanishing of the Jacobian is equivalent to the 
condition for the existence of a singular locus for the 
fanily of streamlines 


(8,4) = Constant 


This singular locus consists of points at which the 
strea -mlines undergo an abrupt change of curvature and 


1¢ans, physically, that the acceleration q oo. of a 
5 


me 
fluid particle is infinite at such points. 


Both Ringleb and Tollmien have shown that the 
singular locus for the streamlines is also the envelope 
of the Mach lines in the plane of flow. The Mach lines 
are related to the streamlines in such a way that the 
component of the fluid velocity normal to a Mach line is 
equal to the local velocity of sound. ‘The Mach lines 
are identical with the so-called characteristic curves 
of the second-order partial differential equations 
for fg and w and are the integral curves of the ordi- 
nary differential equation 
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or 
- 1 aq (56) 


The real solutions of this differential: equation inter- 
preted in the physi¢al’.x, y \ plane.yield, the Mach lines 
for a given flow. The solution of equation (56) is 


fe - 
=o 
8 = 6,22 —— tan oe ie = aisle tan “Life 2 sx (57) 
T 
a J 
where T_ = eae and where 9. issumes the values 


e) 
Py el 


@e 9 along the M= 1 line for a given flow. 


tis recalled oat the Pune taon 
Bo Ss dem.) +, hia) 


introduced in part I (reference 1} in connection with 


the géometric=-mean tyre cof velocity correction formula, 
PSeesoliGLon Oo, the differortial equation 


in the subsonic renge. Observe thet a continuation of 
the function H into the supersonic renge\is given by 
equation (56) as 


evil 7 
qg§ = dq 
g 
In the supersonic range, the function H = 8 - 8, can 


thus be interpreted as ths hodograph of the Mach lines 
for a.given J low. 
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The differential Line elements dx .and dy for the 
Mach lines aes the physical plane are now given. Fron 
equations (} 7) and (56), the line Fieuenta along a Mach 
line for a given flow are 


= 22 (+f 1 cos 9= sin e) (Gus Le <2 


Po (aren 
dy = — (til -1 sin @ +005 6)( 2H tay ms ‘ag 
yok (sus 2lPH By) 


Singular points along a Mach line are characterized by 
the vanishing of the common factor of equations (58) 


~ 


: (58) 


Sheivir -1 = 0 (59) 


OO’ 
4Q 


Equation (59) represents in the plane of flow two pos- 
sible singular loci or “limit lines" for the two families 
of Mach lines associated with the plus end minus signs 
in equation (56). Clearly, the two singular loci cannot 
occur simultaneously since the two conditions cannot be 
satisfied simultaneously. Observe that equation (59) is 
equivalent to the vanishing of the Jacobian given in 
equations (55). Thus, the vanishing of the Jacobian is 
not only the condition for the existence of a singular 
(cusp) locus for the streamlines but also the condition 
for the existence of a limit line (envelope) for the 
Mach lines. 


The existence of a singular locus may be looked upon 
as being equivelent to the vanishing along a curve of 
the Jacobian J oe, of the transformation from the 

> 

hodograph-plane veriables @ and q _ to the physical- 
plane variebles x and y. It-is remarked that 
singular solutions exist for which the Jacobian s(&2 3) 
of the transformetion from the physical-plane veriables x 
and y to the hodogreph-plane variables 8 and q 


vanishes identically in a region of the physical plane. 
In this case, as Tollmien pointed out, 8 and q are 
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no longer independent variables and the flow cannot be 
described in the hcdograph plane. Examples of these 
“missed flows" are the solutions of Meyer (reference 11) 
for supersonic flow inside and outside sharp angles. 


he LS of Special interest to apply the condition 
for the vanishing of the Jacobian to the particular 
solutions pi. and Wi treated in the early part of 
this paper. The expression for the Jacobian for a 
perticular solution 


y 


Po cos ke 
K a 


Wy = =Q° sin k6 
k ite 


is, with the use of equations (5), for k #0, 


eke 
eS e. gin ke + (22) (2-1?) hig costa (60) 
Q- LL Si 


Clearly, this expression for J is positive in the sub- 
Sonic range M<i1. At the sonic value M=1, P, 0 
(see table following equation (23)) and J is again 
positive. At the first zero of P, in the supersonic 
range M>1, Q, # Os hence, "J sis’ negative The 
values of M, for all the pairs of values 8, M for 
which the Jacobien J vanishes, therefore lie between 
M=1 andthe value of M at the first zero of FP, 


(or S;) in the supersonic range. 


By means of the relation 


hd 


e) 
ee a, 
a lhe ta 
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ine) 


if 
the vanishing of the Jacobian yields 


5 
cot ke =F Ls (61) 


ir - 1 


Equation (61) is the relation for pairs of values 9, M, 
which interpreted in the physical =x, y lene constitute 
the limit line for the particular flow Kn? Veo The 
values of M that setisfy equation (61) accordingly lie 
between M= 1 and the value of M at the first zero 

oi Si in the supersonic range. 


This paper is closed with the following remarks on 
limiting values of M in connection with the use of 
velocity correction formulas. The limiting local values 
of M in the case of uniform flow past a prescribed 
boundary, in general, depend on shape parameters. The 
use of a velocity correction formula, however, yields a 
constant limiting value oF M that depends only on the 
particular correction formula used. The geometric-mean 
correction formula yields the value M= 1; the approxi- 
mation of Temple and Yearwood yields M = 1.455; and the 
arithmetic=-mean correction formula given in part I 
(reference 1), which is based on a linear combination of 
a source (limiting value M=1) and a vortex (limiting 
value M=©) or a spiral flow, yields the yalue M = 1.15. 
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